Age Determination Method of Pre-Main Sequence Stars 
with High-Resolution /-Band Spectroscopy*^ 



Yuhei Takagi, Yoichi Itoh 
Graduate School of Science, Kobe University, 1-1 Rokkodai, Nada, Kobe, Hyogo 657-8501 

takagi@stu.kobe-u.ac.jp 
and 
Yumiko Oasa 

Faculty of Education, Saitama University, 255 Shimo-Okubo, Sakura, Saitama, Saitama 338-8570 

(Received ; accepted ) 
Abstract 

We present a new method for determining the age of late-K type pre-main se- 
quence (PMS) stars by deriving their surface gravity from high-resolution /-band 
spectroscopy. Since PMS stars contract as they evolve, age can be estimated from 
surface gravity. We used the equivalent width ratio (EWR) of nearby absorption lines 
to create a surface gravity diagnostic of PMS stars that is free of uncertainties due to 
veiling. The ratios of Fe (8186.7A and 8204.9A) and Na (8183.3A and 8194.8A) ab- 
sorption lines were calculated for giants, main-sequence stars, and weak-line T Tauri 
stars. Effective temperatures were nearly equal across the sample. The Fe to Na 
EWR (Fe/Na) decreases significantly with increasing surface gravity, denoting that 
Fe/Na is a desirable diagnostic for deriving the surface gravity of pre-main sequence 
stars. The surface gravity of PMS stars with 0.8 M is able to be determined with 
an accuracy of 0.1-0.2, which conducts the age of PMS stars within a factor of 1.5, in 
average. 
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1. Introduction 

The age of pre-main sequence (PMS) stars is a vital parameter for studying the evolution 
of young stars. Accurate ages are needed to describe the evolutions of both the photosphere 
(e.g., rotational processes, chemical evolution) and the circumstellar material (e.g., circumstel- 

* Based in part on data collected at Subaru Telescope, which is operated by the National Astronomical 
Observatory of Japan. 

t This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. 
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lar disk, protoplanets, outflows). Determining the age of PMS stars lets us study the evolution 
of the stars themselves and the history of star formation in their star-forming region (e.g., 
Palla & Stahler 2002). Age is usually determined by comparing the luminosity and the ef- 
fective temperature (T e ff) estimated from broadband photometry or spectroscopy with stellar 
evolution tracks on the Hertzsprung-Russell (H-R) diagram (Cohen & Kuhi 1979; Strom et 
al. 1989; Kenyon & Hartmann 1995). However, a precise estimate of PMS star luminosity is 
difficult. This difficulty arises from extinction, uncertainties in the distance measurements, and 
continuum excess caused from accretion and the heated circumstellar disk (e.g., veiling). These 
uncertainties may lead to an incorrect luminosity value. 

The age of a PMS star is able to be determined by estimating the surface gravity (g) 
instead of the luminosity. The surface gravity of a PMS star increases with time because 
of photospheric contraction. Spectroscopic observations are suitable for determining stellar 
parameters such as T e s, log g, metal abundance, and radial velocity, since the absorption line 
profile depends on these parameters. Although the spectra of PMS stars are veiled by the excess, 
the equivalent width ratio (EWR) of absorption lines is less contaminated by veiling (Meyer et 
al. 1998). Meyer et al. (1998) identified the EWR of OH (1.6892/mi) and Mg (1.5760/im) as a 
diagnostic of T cG . They also used the EWR of CO (1.6207/xm + 1.6617/xm) and Mg (1.5760/im) 
as a luminosity indicator (which reflects the surface gravity). Doppmann & Jaffe (2003) derived 
the T e s, radial velocity, and the amount of veiling by fitting observational 2.2\xm spectra with 
synthesized spectra. They estimated log g by calculating the EWR of the Na interval at 2.2/im 
(including the absorption lines of Na, Sc, Si) and the CO bandhead at 2.3/im. When the EWR 
is calculated by taking the ratio of absorption lines separated, a correction for the amount of 
veiling is necessary, as it is wavelength dependent. 

To derive the EWR directly from observed spectra without any veiling correction, we 
calculate EWRs using closely spaced absorption lines. This is practicable because the effects 
of veiling are nearly constant over small wavelength ranges. We present a new method for 
determining the age of mid- to late-K type PMS stars using EWRs of nearby absorption lines 
obtained from high-resolution /-band spectroscopy. We observed 30 giant stars, 6 main-sequence 
stars, and 4 weak-line T Tauri stars (WTTSs) with T e g of around 4200K to derive a relationship 
between log g and the EWRs of Fe (8186.71, 8204.91) and Na (8183.31, 8194.81). 

2. Method 

To make highly accurate determinations of PMS star age, precise log g estimation is nec- 
essary. We selected the wavelength range in /-band (8180! to 8210A), where strong atomic lines 
exist. This wavelength is suitable for low-mass PMS stars of which T e g is around 4000K. Such 
PMS stars in the Hayashi phase (Hayashi 1961) evolve with roughly constant T e g. Therefore, 
their age may be approximately determined from log g. 

The absorption lines used in the EWR must be selected with care to minimize uncertainty 
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in the resultant log g value. We selected the J-band Na lines (8183. 3A, 8194. 8 A) and Fe lines 
(8186. 7 A (blended with vanadium), 8204. 9A). These lines are sensitive to both T e g and log g. 
Na lines are one of the strongest lines in the /-band. They are reached to the limiting depth, and 
the wing of these absorption lines is pressure-broadened in stars with large g. Fe lines are also 
strong lines in the J-band, but they are not reached to the limiting depth. The Fe line strength 
decreases with increasing g, depending on the growth of continuum opacity. Therefore, EWRs 
calculated from Na and Fe absorption line pair (Fe/Na) are considered suitable for deriving 
stellar surface gravities. In addition, assuming the abundance ratios of the elements are equal, 
EWR has no dependence on metal abundance. 

However, the relation between Fe/Na and surface gravity is yet unknown. To derive this 
relationship, a high-resolution spectroscopic observation of stars with fixed T e g and varying 
log g was necessary. Since the surface gravities of PMS stars are between that of giants and 
main-sequence stars, we observed giants, main-sequence stars. Also the WTTSs were observed 
to verify the effectiveness of this method to PMS stars. 

3. Observations and Data Reduction 

3.1. Object Properties 

The selected sample objects and their properties are listed in table 1. We selected 30 
giants (including subgiants) and 6 main-sequence stars to estimate the relationship between g 
and the Fe/Na. All targets are single stars. 25 giants and 4 main-sequence stars have T e s of 
4200-4300K, since the radiation peak of these stars is in the /-band and they are bright enough 
to observe. Five giants and two dwarfs with T c g of 4600-4700K were included to provide an 
estimate of the dependence of EWRs on T e e. 

To derive the surface gravity in advance, giants and main-sequence stars of which parallax 

was measured by the Hipparcos satellite (Perryman et al. 1997) were selected. The log g of 

these objects was calculated as follows: 

M T L 

log# = log— + 41og— -log— + log# , (1) 

1V1 & 1 L/q 

where M is the mass, T the effective temperature, and L the luminosity. The luminosity was 
calculated from the l/-band magnitude given in SIMBAD (taken from previous photometric 
studies), and the parallax. For the effective temperature, the T e ff-spectral type relationship 
(Lang 1992) was used. The error of the effective temperature is not included in the log g esti- 
mation and therefore the error of the log g may be underestimated. The relationship of M and 
spectral type in Drilling & Landolt (2000) was used to derive the masses of the main-sequence 
stars. The giant star masses were estimated by comparing luminosities and temperatures with 
the evolutionary tracks given in Lejeune & Schaerer (2001) (figure 1). Since giants smaller than 
2.0 M & evolve rapidly after the helium flash, we used the pre-flash evolutionary track for these 
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Fig. 1. The H-R diagram for the giants and main-sequence stars. The dashed lines are the evolutionary 
tracks from Lejeunc & Schaerer (2001); M = 0.90, 1.00, 1.25, 1.50, 1.70, 2.00, 2.50, 3.00, 4.00, and 5.00 
Mq from the bottom. The solid line indicates the isochrone for 10 9 yr (Baraffe et al. 1998). 

stars. The derived log g of /3 Gem (2.73±0.01) is consistent with the result of a spectroscopic 
study (Drake and Smith 1991), in which it was estimated as 2.75±0.15. The large error in the 
surface gravity of the giants is attributable to the errors in luminosity (figure 1), which arise 
from the parallax errors. 

Four WTTSs belonging to the Taurus-Auriga star-forming region were also added to 
our samples to test the effectiveness of the derived EWR-log g relationship. The log g of 
these WTTSs were also estimated from equation (1). Temperatures were estimated using Lang 
(1992), following the procedure for the main-sequence stars. We derived luminosities from 
photometry, including corrections for extinction (Ay), since the excess from the circumstellar 
disk is small. The Ay for HBC 374 and V827 Tau were taken from Kenyon & Hartmann 
(1995). For the other objects, RX J0452.5+1730 and RX J0459.7+1430, we estimated the A v 
by fitting the SED corrected by an arbitrary Ay to the V, J, H, and K-band magnitudes. 
The V^-band magnitude was taken from SIMBAD, and the 2MASS catalog was used for the 
J, H, and i^-band magnitudes. The estimated Ay was 0.69 and 0.67 for RX J0452. 5+1730 
and RX J0459.7+1430, respectively The mass of each star was estimated by using the H-R 
diagram, comparing the temperature and luminosity with the evolutionary model of Baraffe et 
al. (1998). The mass of HBC 374 and V827 Tau were 0.8 M & . On the other hand, the masses 
of RX J0452.5+1730 and RX J0459. 7+1430 were 1.20 M and 1.25 M , respectively. The 
distance to the Taurus- Auriga molecular cloud was fixed at 142±14pc (Wichmann et al. 1998). 
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3.2. Observations 

High-resolution /-band spectroscopy of 30 giants and 6 main- sequence stars was carried 
out on 2007 January 17 to 24 and 2008 March 24 to 27 using the High Dispersion Echelle 
Spectrograph (HIDES) on the Okayama Astrophysical Observatory 1.88 m telescope. The slit 
width was set to 220/mi, giving a resolution of ~60000. The angle of the red cross-disperser was 
selected to center the 69th order of the echelle spectrograph (8 165^8300 A) on the 2048x4096 
CCD. The integration time range was 90^7200 sec, yielding a S/N of ~100. The maximum 
integration time was set to 1800 sec to reduce the risk of cosmic rays. Exposures of faint objects 
were divided into several frames. Flat frame were taken using the flat lamp equipped in the 
instrument, and were obtained at the beginning and end of each night's observations. Bias 
frames and Th-Ar lamp frames for wavelength calibration were taken with the same frequency. 

The spectra of the 4 WTTSs were obtained using the Subaru Telescope with High 
Dispersion Spectrograph (HDS) on 2007 September 18. The slit width was set to 0.6" 
(R~60000). The instrument was set to the StdNIRb mode, to obtain the Na and Fe ab- 
sorption lines in the 73rd order. The integration time for each object was 900 sec. Flat, bias, 
and Th-Ar frames were also obtained. 

3.3. Data Reduction 

The data were reduced using Image Reduction and Analysis Facility (IRAF) software 
package 1 . The same procedures were applied to both HIDES data and HDS data. First, the 
bias was removed for each frame by subtracting the average level of the overscan area of each 
frame, and also by subtracting the averaged bias frame. Flat fielding was accomplished with 
the normalized flat frame. Cosmic rays and scattered light were removed. Next, the spectrum 
was extracted from each object frame. The comparison data obtained with the Th-Ar lamp 
were used for wavelength calibration. When the object was taken several times, the spectra 
were combined after wavelength calibration to improve the signal-to-noise ratio. The spectra 
were normalized at the continuum level in the last procedure. Spectra samples are displayed 
in figure 2. 

The SPLOT task was used to measure the equivalent width (EW) of each line, and the 
line was fit with the Voigt function. Since telluric lines appear in this wavelength range, we 
de-blended these lines from the photospheric lines. The EWs of some Na and Fe lines were 
difficult to measure because of heavy telluric line blending. The error in EW was estimated 
from uncertainty of the continuum level. The measured EWs of each object are listed in table 
2. 

1 IRAF is distributed by the National Optical Astronomy Observatory. 
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Fig. 2. Sample spectra of a giant, a WTTS, and a main-sequence star. The Doppler shift of each object 
is corrected. The absorption lines used to derive the EWRs arc marked. Other sharp lines are mostly 
telluric lines. The Na (8183.31) and Fe (8186.7A) in the WTTS spectrum are blended with the telluric 
lines. 

4. Results and Discussion 



4-1. Relationship between EWR and Surface Gravity 

We derived a relationship between log g and the EWR of Fe and Na absorption lines 
(Fe/Na) for 25 giants and the 4 main-sequence stars of which T e g are around 4200-4300K. We 
calculated the EWR using four line pairs: Fe (8186.71) / Na (8183.31), Fe (8186.71) / Na 
(8194.81), Fe (8204.91) / Na (8183.31), and Fe (8204.91) / Na (8194.81). Results are shown 
in figure 3. Filled squares and triangles mark the EWRs of the giants and main-sequence stars, 
respectively. The solid line in each graph is the best-fit line for these objects. In every case, 
the line ratios decrease significantly with increasing log g. The averaged slope of the curves 
in log g = 3.5-4.0 (a typical surface gravity for PMS stars) are -0.085, -0.067, -0.084, -0.064 
in Fe (8186.71) / Na (8183.31), Fe (8186.71) / Na (8194.81), Fe (8204.91) / Na (8183.31), 
and Fe (8204.91) / Na (8194.81), respectively. These EWRs are considered to be an excellent 
diagnostics for deriving the surface gravity. 

Next, we calculated the EWRs of five giants and two main-sequence stars of which T e g- 
is 4600-4700K. The EWRs of these objects are shown as open squares and triangles in figure 3. 
The best-fit line in each EWR is indicated as a dashed line. The EWRs increase for higher- Teg- 
objects at constant gravity over the range of log g = 3.5-4.0. A comparison of fitted curves 
indicates that the EWRs of high-T c g- objects are 0.007-0.032 larger than the low-Teg EWRs in 
log g = 3.5. 

6 




012345 012345 
log g log g 



Fig. 3. EWR versus log g (g: cm/s 2 ). The squares mark the EWRs for the giants. The triangles mark the 
EWRs for main-sequence stars. The filled plots denote the EWRs of objects with T e g around 4200-4300K, 
and the open plots are for 4600-4700K objects. The solid lines are the fitted curve for the 4200-4300K 
objects. The lines for the 4600-4700K objects are shown as dashed. 
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4-2. Model Atmospheric Spectra Comparisons 

We synthesized spectra using the SPTOOL program (Takeda 1995), which calculates 
model spectra based on the ATLAS 9 model (Kurucz 1993). The model spectra were calculated 
for T e ff of 4200K and 4750K, and a range in log g of 1.0-5.0, and the EWRs were estimated 
(figure 4). All of the EWRs calculated from the synthesized spectra decreased with increasing 
log g. This is consistent with the EWRs for the observed spectra. 

In addition, although the EWRs of the observed spectra and the data from the SPTOOL 
are poorly matched especially in the EWRs from the 4750K model spectra and the Fe (8204. 9 A) 
ratios, the EWRs of the high-T e g objects are larger than those of low-T e g objects with a log g 
range larger than 3. This occurs because the EW of the Fe lines increases, and the wing of the 
Na lines decreases in high-T e ff spectra. 

This discrepancy in observed spectra and model spectra could arise from improper values 
of the excitation potential and the damping constant of the Fe and Na lines. Since the both 
observed and synthetic EWRs match in the Fe (8186.71) / Na (8183.31) and Fe (8186.71) / Na 
(8194.81) EWRs, the excitation potential and the damping constant of Fe (8186.71) line and 
both Na lines are seem to be estimated correctly. In the Fe (8204.91) case, since the observed 
EWs were weaker than the modeled EWs and has no strong wings, the excitation potential 
may be underestimated. Further improvements on the line data will increase the utility of this 
log g diagnostic. 

4.3. The EWR of WTTS 

We next compared the EWRs of WTTSs with the EWRs of giants and main-sequence 
stars. Since the Na (8183.31) and Fe (8186.71) lines of the WTTSs were blended with the 
telluric lines, only the Fe(8204.9l) / Na (8194.81) was obtained. 

As seen in figure 5, the EWRs of two WTTSs (RX J0452.5+1730 and RX J0459.7+1430, 
marked with open circles) with T c q =4590K are larger than the EWR of the WTTS (HBC 374, 
marked with filled circle) of which T e g is 4060K. This is consistent with results from the fitted 
curves and the model spectra mentioned in section 4.2. The EWRs of high-Teg WTTSs agree 
closely with the curve from the 4600-4700K objects, and the EWRs of HBC 374 and V827 Tau 
fits the curve approximated by the 4200-4300K objects. This result indicates that the surface 
gravity of PMS stars is able to be determined from this relationship. 

By using this relationship, the log g of these two WTTSs are estimated at 3.95, which 
correspond to the age of 2.8 Myr in the evolution model of Baraffe et al. (1998). This result 
agrees with the age derived from the photometric observation using the extinction amount of 
Kenyon & Hartmann (1995) (table 3). However, using the amount of extinction of Furlan et 
al. (2006) results the log g of HBC 374 and V827 Tau to 3.263 and 3.627, respectively. The 
ages of these two WTTSs estimated by using this Ay will be less than 1 Myr, which mismatch 
our result. 
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Fig. 4. EWR versus log g (g: cm/s 2 ) with the EWR of the model spectra created by SPTOOL. The 
squares indicate the EWR of the giants, and the triangles represent the main-sequence stars. The filled 
plots denote the EWRs of the 4200-4300K objects, and the open plots are for 4600-4700K objects. The 
plus and cross signs are the EWR from model spectra of T e g = 4750K and T e s = 4250K, respectively. 
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Fig. 5. EWR (Fe (8204.91) / Na (8194.81)) versus log g (g: cm/s 2 ) for the EWR of WTTSs. The 
squares, triangles, and the fitted curves are from figure 3. The open and filled circles indicate the EWRs 
of WTTS of which T cS = 4590K and T cff = 4060K, respectively. 
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4-4- Age Determination Accuracy 

The age of a PMS star is able to be estimated by comparing the derived surface gravity 
with an evolution model (Baraffe et al. 1998; Siess et al. 2000). In our method, the accuracy 
of the calculated ages depends on the precision of the surface gravity value obtained for the 
PMS stars. This is dependent on the error in the EWR. We assumed the EWR error of the 
PMS star as 7.0 x 10~ 3 , which nearly corresponds to the error of HBC 374. The log g error 
calculated from this uncertainty is 0.1-0.2 in the 3.5-4.0 range. Consequently, by comparing 
with an evolution model of a 0.8 M star in Baraffe et al. (1998), the age of the PMS star can 
be determined within a factor of 1.5, in average. This precision will be improved by deriving the 
general EWR error in the other three EWRs. To adapt this method to PMS stars of different 
temperatures, accurate estimations of the relation of EWR-log g relationship for each T e g range 
are required. In addition, an independent T e g estimation is necessary in order to determine the 
log g of the pre-main sequence star using this technique. 

5. Conclusion 

We carried out high-resolution optical spectroscopy to create a new age determination 
method for late-K type PMS stars. To obtain a surface gravity indicator from the spectrum, 
we derived the EWR using nearby absorption lines, which avoid problems from veiling contam- 
ination and distance uncertainty. 

The EW of the Na (8183.3A and 8194.8A) and Fe (8186.7A and 8204.9A) absorption lines 
in the /-band was used in this work. To derive the relationship between Fe/Na and surface 
gravity, we observed 25 giant stars, 4 main sequence stars, and WTTSs with the Okayama 
Astrophysical Observatory 1.88 m telescope and HIDES, and the Subaru Telescope and HDS. 
The T c g of the objects were 4200-4300K and the surface gravity was known from the previous 
photometric studies. All four EWRs decreased with increasing surface gravity. Therefore, they 
are efficient diagnostics of surface gravity. Using the Fe (8204. 9A) / Na (8194. 8A) EWR, the log 
g of PMS stars with 0.8 M can be estimated with an uncertainty of 0.1-0.2. From comparisons 
with an evolution model, we concluded that their ages can be determined within a factor of 
1.5. 

Five giants, two main sequence stars, and two WTTSs with higher- T c g (4600-4700K) 
were also observed to estimate the temperature dependence of the EWR. The EWRs of the 
high-Teg objects were larger than those of the low T cS in the typical PMS star surface gravity 
range. Therefore, to accurately determine the ages of PMS stars with masses larger or smaller 
than 0.8 M Q , new EWR-log g relationship for each T e g range are required. A precise estimation 
of EWR-log g relationship allow us to determine the accurate age of PMS stars and to study 
the evolution process of star and circumstellar materials. 



11 



We are grateful for the assistance given by the staff of the Okayama Astrophysical 
Observatory during the observations. We also thank the staff of the Subaru Telescope for the 
service observation. 

References 

Baraffe, I, Chabrier, G., Allard, F. & Hauschildt, P. H. 1998, A&A, 337, 403 
Cohen, M. & Kuhi, L. V. 1979, ApJS, 41, 743 
Drake, J. J. & Smith, G. 1991, MNRAS, 250, 89 
Doppmann, G. W. & Jaffe, D.T. 2003, AJ, 126, 3030 

Drilling, J. S. & Landolt, A. U. 2000, in Allen's Astrophysical Quantities, ed. A. N. Cox (New York: 

Springer- Verlag) , 389 
Furlan, E., et al. 2006, ApJS, 165, 568 
Kenyon, S. J. k Hartmann, L. 1995, ApJS, 101, 117 
Hayashi, C. 1961, PASJ, 13, 450 

Kurucz, R. L. 1993, Kurucz CD-ROM, No. 13 (Cambridge: Smithsonian Astrophysical Observatory) 
Lang, K. R. 1992, Astrophysical Data: Planets and Stars (New York: Springer- Verlag) , ch. 9 
Lejeune, T. k Schaerer, D. 2001, A&A, 366, 538 

Meyer, M. R., Edwards, S., Hinkle, K. H., k Strom, S. E. 1998, ApJ, 508, 397 

Palla, F. k Stahler, S. W. 2002, ApJ, 581, 1194 

Perryman, M. A. C, et al. 1997, A&A, 323, L49 

Siess, L., Dufour, E. k Forestini, M. 2000, A&A, 358, 593 

Strom, K. M., Strom, S. E., Edwards, S., Cabrit, S. & Skrutskie, M. F. 1989, AJ, 97, 1451 
Takeda, Y. 1995, PASJ, 47, 287 

Wichmann, R., Bastian, U., Krautter, J., Jankovics, I. & Rucihski, S. M. 1998, MNRAS, 301, L39 



12 



Table 1. The list of the observed stars and their properties. 



Name 


Sp. Type 


Teff (#) 


log 9 {g- cm s 2 ) 


S/N 


Giants 


p Gem 


KOlllb 


4750 


O 70 1 A A1 

2.73±0.01 


1 CA 

160 


TTT~\ 1 A/ 1 r 

HD 19656 


KOIII 


4750 


O 1 1 A AO 

2.38±0.02 


100 


HD 7oz91 


K11V 


4600 


o 70+0.02 

A - 1 z -0.04 


1 TA 

170 


HD l(J(Jz(J4 


K11V 


A O AA 

4600 


o m +0.19 


1 /i A 

141) 


HD 145148 


K 1.51V 


4530 


O OO 1 A AO 

3.23±0.03 


1 TA 

170 


HD 45512 


K2111-1V 


4420 


O f r* 1 A AO 

2.55±0.03 


130 


HD 101978 


K2IV 


4420 


i co+0.27 

J-^-i.o? 


160 


TT"T\ -1 A/11 AO 

HD 106102 


K2111-1V 


4420 


o ci +0.08 

2-61_ .o5 


140 


TTT~\ 1 AO OAA 

HD l(Joz99 


K21V 


A A OA 

4420 


O /IG+0.11 

2 - 46 -0.13 


1 AA 

100 


T TTA -1 I AC A1 

HD llUoUl 


TVOTTT TT J 

K2111-1V 


4420 


O A 1 I A A A 

z.91±l).(J4 


1 1 A 

110 


TTT~\ 1 fl1 1 AC 

HD 121146 


K2IV 


4420 


o r A I A A o 

2.50±0.02 


190 


TT"T\ -1 A >-7 1 /JO 

HD 147142 


K2IV 


4420 


o ao+0.06 

2 -03I .o3 


110 


HD 49738 


K3III 


4200 


i qo+0.20 
l.OO_ 25 


110 


HD 65759 


T/OTTT 

K3111 


/I OAA 

4200 


1 GG + 0.09 

1.66+ n 


1 1 A 

110 


HD 86369 


K3III 


4200 


1 on+0.06 


160 


TTT"\ OOOOI 

HD 88231 


K3III 


4200 


1 *7C + 0.10 


80 


TT"T\ -i AOOOO 

HD 102328 


K3III 


4200 


O 1 r 1 A AO 

2.15±0.02 


90 


TTT~\ 1 AAA1 O 

HD 109012 


K3111-1V 


4200 


o oc+0.12 
AZ °-0.14 


160 


TTT\ 1 i op on 

HD lldbol 


K3111 


A OAA 

4200 


1 CG+0.11 

i-^-o.is 


1 OA 

120 


TTTA 1 1 r /I TO 

HD 115478 


K3III 


4200 


O A1 I A AO 

2.01±0.03 


140 


HD 118839 


K3III 


4200 


1 O O I A 1 A 

1.83±0.10 


130 


HD 129245 


K3III 


4200 


o ao+0.03 

2.02+0.02 


90 


I I n 1 QOOA/i 
11 U 10Z0U4 


ivolll 


/i°nn 

4ZUU 


i 7O+0.16 
x - ' °-0.13 


i qn 

loU 


HD 137688 


T/OTTT 

K3111 


A OAA 

4200 


1 7A+0.17 


1 OA 

120 


HD 146388 


K3III 


4200 


z-UD_o.o4 


140 


HD 156093 


K3III 


4200 


-, ? r+0.10 

L - ' °-0.06 


170 


HD 163547 


K3III 


4200 


i OO+0.05 
± -°°-0.06 


190 


HD 10824 


K4III 


4000 


i ro+0.08 
^"''-O.OQ 


150 


HD 20644 


K4III 


4000 


-i 97+0.04 
-"-•^ ' -0.08 


100 


HD 23413 


K4III 


4000 




120 


Main-sequence stars 


GL 105 


K3V 


4730 


4.46±0.01 


270 


GL 183 


K3V 


4730 


4.44±0.01 


110 


GL 141 


K5V 


4350 


4.37±0.01 


120 
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Table 1. (Continued.) 



Name 


Sp. Type 


T off (K) 


log 9 cm s" 2 ) 


S/N 


GL 380 


K5V 


4350 


4.85±0.00 


220 


GL 397 


K5 


4350 


4.73±0.02 


130 


GJ 3678 


K5V 


4350 


4.30±0.02 


150 


WTTS 


RX J0452.5+1730 


K4 


4590 


a oo +0.05 
^■ zz -0.07 


60 


RX J0459.7+1430 


K4 


4590 


4 TQ+0.07 
^• io -0.08 


60 


HBC 374 


K7 


4060 


o no+0-08 
°- JZ -0.09 


70 


V827 Tau 


K7 


4060 


J - Jz -0.10 


70 
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Table 2. The list of the measured EW for each object. 





Name 




EW (1) 








Na (8183.31) 


Fe (8186.71) 


Na (8194.81) 


Fe (8204.91) 


Giants 


P Gem 


n oan+0 009 


A 1AA_LA A AC 

L).1(J(J±L).(J(Jd 


n 07/1 +0.008 
0-374_ .oo9 


A AO/I _I_A AA1 


HD 


19656 


n OR/1+ 014 
0.264_ 0010 




A Q07+O.OI6 
U.OO ( _ .007 


n in/1 i n nnQ 


nn 
rlU 


7R001 

/ ozy i 


0.251±0.006 


« ^ _ . 4_n 007 


0.342±0.008 


rAri1 4-0 003 

0.091^;^ 


HD 


1 ( i ( 1 .''1 ( 1 A 

100204 


n 0A 7 +0.007 


n -1 1 .,+0.007 
U.114_ 008 


A OCQ+0.007 

u.ooy_ 009 


U. lU4ztU.UUz 


HD 


145148 


n oi O+0.006 

0-313I 005 


0.072±0.00l 


0.43l±0.006 


0.069±0.002 


HD 


a tr c 1 o 

4551z 


n on /i+0 014 


nil o+O 006 
0-H3Io.005 


0.434±0.008 


0.111±0.003 


111*1 


iuiy ( o 


A ota+0.006 

U.Z(»_g QQg 


n 1 ec +0.006 

0-165I .oo5 


n 0/1-7+O.OI2 

0.347+0.013 


A 1 /i A I A A A O 

0.149±0.008 


HD 


1 A A 1 AO 

lUolUz 


n q/fi+o.oii 
0.341_ 006 


n 1 qn+0001 

U. 10U_q 


A /IQR+0- 017 
0.49b_ .oi9 


A lAQ+O- 001 

u.iuy_ 002 


HD 


1 AOOAA 

108299 


n O7/1+0.007 


A 1 OO I A A A 1 

0.138±0.001 


n oon+0.007 
u.ozy_ 006 


n 1 on+0.004 

0-130I .oo5 


HD 


110501 


A 01 C I A AAA 

0.315±0.006 


n 1 nc+0.005 

0-105I .oo9 


n /(cn+0.013 

0.460+0.012 


A AO/ 1 1 A AAO 

0.086±0.00z 




i o i 1 a a 

12114b 




t\ -1 rtr»4-0 006 

0-l26±£;5o5 


„ 4-0 002 
u ' oyo -0.001 


r , 1 ^rt +0 001 

0.103t°;X 


HD 


1/171/10 
14 ( 14z 


n oqq+0.006 
u.zyo_g 


A 1AA 1 A A AC 

0.109±0.005 


A OHO 1 A AAO 

0.393±0.008 


n 1 99+0.000 
u.izz_ 001 


HD 


A n70 O 

49 ( oo 


u.ouy±u.uii 


n 1 90+0. 007 
u.izy_ 005 


A O.7C+0.004 
u -°'°-0.005 


U.10U±U.UU1 


HD 


65759 


n oqq+0.009 
u.zyo_ 008 


n 1 n nn 1 

0.114±0.UU1 


n ocn+0.010 

0.360I .oo7 


n 1 oc+0.000 

0.136I 00 i 


iii^ 

HD 


ooooy 




n 1 /i e+0.006 
°- 145 -0.005 


A ^Sl +0 007 
U.oSl_o 008 




HD 


oozol 


n oni _i_n m o 
L).291±U.013 


A 1 O _1_ A A A /I 

U.lzo±U.(J(J4 


n qcc+0-014 
0.368_ .oi5 


U.131±U.001 


HD 


1 ao o o o 

102328 


n Q7/I+0.005 

0.374_ 0012 


n i c;o.-i-n nno. 
U.IOUitU.UUo 


U.<34z±U.Ulz 


n 1 O/l +0.002 

0.134_ 00 o3 


HD 


109012 


A OOA I A AA/ 1 

0.289±0.00o 




A 900+0.010 

u.ooo_o 012 


A 1 AO I A AAO 

0.108±0.003 


HD 


1 1 0/J07 


rv r»o^+0 003 

0.286™2 


~ -4-0 008 
0-H5-0.004 


^ ^j_q-)-0.004 
U.008_o.oo7 


A , , o+O 004 

0-142±o^o6 


TTT^ 

rlU 


1 104 ( o 


n on/i+0.002 

0.294r o . oi 


n 1 o/i +0.003 
0-124-0.002 


A /] O9+0.004 
U. z lOZ_ 002 


A 1 9C+0.004 
U.1ZO_ 005 


HD 


118839 






A Q77+0.009 
u -° ' ' -0.010 


A 1 1 04-A AAO 
U.llzitU.UUZ 




i ono /i £ 
lzyz4o 






0.398±0.006 


0.127±0.004 


HD 


132304 




n 1 QO_i_n nno 
U.loz±U.UUz 


a onA+0 003 
0.300_o. 00 4 


n 1 on_i_n nno 
U.120±U.00o 


HD 


137688 


0.248™ 


n 11O+0.007 
U.llO_ 003 


A Qrn+0.013 

u.oou_ 010 


n i 24 +0.002 
U.1Z4_ 004 


HD 


146388 


n 977+0.004 
u - z ' ' -0.005 


0.110±0.007 


0.394±0.001 


0.106±0.001 


HD 


156093 


263+ 0012 


n i 42 +0.003 

u.i z ±z_ 002 


0.343±0.001 


0.131±0.001 


HD 


163547 


2 64+ 006 


n 1 99+O.OO7 
u.izz_ 005 




n 194 + 0.004 

U.1Z1_ 002 


HD 


10824 


S20+ 009 


a 174+0001 

Uml ' ^-0.002 


0.398±0.011 


U.±U( _o.005 


HD 


20644 


n q.47+0.016 

U.O1/_ 014 


0.164±0.004 


407+ - 016 

_ .007 


0.166±0.005 


HD 


23413 


n 099+0.010 
u.ozz_ 006 






n ion+0.009 

u.ioy_ 006 



Main-sequence stars 


GL 105 
GL 183 


0A78t° fd 
0-629™ 


A7A+0 00 5 

U.U(U_ 004 


0.596±0.004 0.046±0.001 
a 065+ 007 
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Table 2. (Continued.) 



Name 




EW (1) 






IN a ^oloo.o/i ) 


Fe (8186.71) 


Na (8194.81) 


rC (ozU4.y/i / ) 


GL 141 


0.82618:855 


1 T 1 + 007 

U.lll_g QQg 


i 171 +0.020 
1 - L 1 i -0.023 


n n 7 o+0.001 
u - u '°-0.002 


GL 380 




U - UJ ' -0.004 


1 Q 4 t;+0.048 
l.O^O_ 052 


066+ 004 
u.uoo_ 002 


GL 397 


o.94ol8;81I 


0.117+0.011 


1 OQ7+0.035 
1 _o.028 


0.069+0.002 


GJ 3678 


0.50618:88? 


n58+ 008 

U.UOO_ .006 


u - ' J'-o.ooe 


036+ 002 

U.UOD_ 0Q1 


WTTS 


RX J0452. 5+1730 






fl 779+0.026 
u - ' ' z -0.033 


0.08218:8^ 


RX J0459. 7+1430 






7 20+ 033 
u. (^u_ 065 


o.o86l8:8f 7 


HBC 374 






ORQ+0 037 
u.»uo_ 026 


0.085+0.007 


V827 Tau 






i 990+O.O2I 
l.zza_ ,053 


o.io9l8:g? 



The lines heavily blended with telluric lines were unmeasured (expressed " — "). 



Table 3. The age of the WTTSs estimated from EWR and photometric method. 



Name 



Method 



EWR 



HBC 374 
V827 Tau 



2.8 
2.8 



Photometric 



KH95* 



0.76 
0.28 



2.5 
2.5 



* Kcnyon k. Hartmann (1995) 
t Furlan et al. (2006) 



F06^ 



;e (Myr) A v Age (Myr) A v Age (Myr) 



2.40 
1.00 



<1.0 
<1.0 
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